Dimensional Crossover Driven by an Electric Field 
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We study the steady-state dynamics of the Hubbard model driven out-of-equilibrium by a constant 
electric field and coupled to a dissipative heat bath. For very strong field, we find a dimensional 
reduction: the system behaves as an equilibrium Hubbard model in lower dimensions. We derive 
steady-state equations for the dynamical mean-field theory in the presence of dissipation. We discuss 
how the electric field induced dimensional crossover affects the momentum resolved and integrated 
spectral functions, the energy distribution function, as well as the steady current in the non-linear 
regime. 



There is a growing interest in quantum physics out 
of equilibrium motivated by new experimental realiza- 
tions in non-linear transport in devices, heterostructures 
and cold atoms [J-Q. On the theoretical side, new tech- 
niques are being developed in the context of correlated 
electrons driven by electric field, ranging from exact di- 
agonalization and density matrix renormalization group 
techniques for treating one-dimensional systems to ex- 
tensions of the dynamical mean-field theory (DMFT) in 
out-of-equilibrium situations Important problems 

such as the suppression of the Bloch oscillations by the 
electronic interactions fs^] and the dielectric breakdown 
of the Mott insulator j9|] have been studied with these 
methods. 

In this paper, we focus on the non-equilibrium steady- 
state dynamics of a correlated metal driven by a static 
and uniform electric field. We formulate a consistent de- 
scription of the steady-state by including a thermostat, 
which plays a key role to prevent the subsequent cur- 
rent to heat up the system. We show analytically that 
the system undergoes a dimensional crossover: for very 
strong fields, the physics reach thermal equilibrium in 
lower dimensions. We derive steady-state equations for 
the DMFT in this driven scenario and simply solve the 
impurity by the so-called iterated perturbation theory 
(IPT) . We study the effects of the electric field and illus- 
trate the dimensional crossover by computing experimen- 
tally accessible quantities including the spectral function, 
the energy distribution function or the steady current. 

Model. We consider the d-dimensional Hubbard 
model on a hypercubic lattice of lattice spacing a. A 
static and uniform electric field is set along one of the 
axes of the lattice: E = Eu^ with E > 0. The La- 
grangian of the system is given by (we set h = 1) 



'Cs = ^ Cicr [idt - <Pi{t)] Ci<j - [/ ^ Ci^Ci^CilCil 



conj . , 



(1) 



where Cia and Cia are the Grasmann fields represent- 
ing an electron at site i with spin a g {tj4-}- U is 
the on-site Coulombic repulsion between electrons and 
tij = (a/27r)'' Jdk e''''^'J e(k) sets the hopping amplitude 
between two neighboring sites. Because of the periodicity 



in the lattice, all integrals in k-space are computed in the 
Brillouin zone {i.e. between —ir/a and tt/o in each direc- 
tion). The Peierls phase factors aij{t) = q J^' dx •A(i, x) 

are required by the gauged U(l) symmetry associated 
with the conservation of the charge q of the electrons. 
We gather the scalar potential and the vector poten- 
tial A in the gauge field A^^ = {(j)/q,A). 

The electric field acts as a source of energy and we 
provide a heat sink by coupling each site i to an inde- 
pendent thermal bath 13, 11|. These are reservoirs of 



non-interacting electrons in equilibrium (and not affected 
by the electric field) at temperature T and chemical po- 
tential /iQ. The latter controls the electronic occupation 
of the system. We choose a gauge-invariant bi-linear cou- 
pling with the system: £sb = lJ2iai e'®'^*^5i^iCi„ + conj. 
where 7 is a coupling constant, the b's represent the 
bath electrons, / labels their energy levels and Oi{t) = 
J dt' (j)i{t'). A regular dissipation in the whole spectrum 
is achieved if the bandwidth W of the reservoirs is larger 
than any other energy scale. 

The energy scales involved are the strength of the hop- 
ping controlled by eq = e(k = 0)/ V^d, the strength of the 
electronic repulsion U, the temperature T that sets the 
scale of the thermal excitations (we set fee — 1), and the 
strength of the dissipation F = The energy scale 

introduced by the electric field is given by Iqjii'a. 

Steady-state equations. — We analyze the dynamics 
by use of the Schwinger-Keldysh formalism • Corre- 
lations and field being switched on far in the past, the 
properties of the steady state do not depend on the de- 
tails of the initial conditions. The steady-state physics is 
invariant in space and time. In particular, a mean- field 
analysis should not depend on the location and time at 
which the impurity is singled out from the bulk. This 
motivates us to work with a covariant formalism rather 
than imposing a particular gauge which would artificially 
break the symmetry of the problem. Steady-state equa- 
tions are then derived -without having to solve the tran- 
sient dynamics- by imposing space and time translational 
invariance. 

To shorten the notations, we drop the spin indices and, 
while working on the lattice, we adopt continuous space 
notation [13]: space-time coordinates are given by 2:^ = 
(t,x). The Schwinger- Dyson equations in the Keldysh 
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basis read: 



D{x,xi) o^^G^{xi,x') = 5{x ~ x') , 

G''{x,x') = G^(x,xi) o,^S^(xi,:r2) o^^G'^ix'.x^y , 

with D{x,x') EE 5(x - x') [idf - <j){x)] + 6{t - t')t{x - 



x')e' 



S (x, a:') and o^.^ is the convolution prod- 



uct (inner summation on xi). and G^ are respec- 
tively the retarded and the Keldysh components of the 
2 by 2 Green's function matrix: = G++ - G+- 

and G^ = i[G+- + G-+]/2. Both the thermal en- 
vironment and the electronic repulsion contribute to 



the self-energy: S^/-^ = E^''^ -h S^'". In equilib- 
rium, Green's functions (and self-energies) are related 
through the fermionic fluctuation-dissipation theorem 
(FDT): G^{k) = - tanh (^i^) Im G«(fc) (the same re- 
lation holds between and ImE^) where fc'^ = (a;,k) 
are the Fourier conjugates of 



^R/K 



After a Wigner transform [2j], the Schwinger-Dyson 
equations are expressed in terms of the center-of-mass 
coordinates X'^ and the conjugated variables k'^ of the 
relative coordinates a;'': 

D{X;k)-kG^{X;k) = 1, 

G^{X; k) = G^(X; k) * E^(X; k) * G^(X; k)* . 

The star product is the non-commutative Moyal product 

defined by ★ = exp {^\c)ki^dx^ — '(ho^dkl]) [Ul- Left (resp. 
right) arrows indicate that the derivative operators act 
on the left (resp. right) and the Einstein summation 
convention is used. 

Schwinger-Dyson equations can be rendered explic- 
itly covariant using the covariant derivatives id^ — 



qAf^{X) and their Wigner transforms 



'{k,X) = 

(zJjK) [15|. Once Green's func- 



tions are re-expressed in terms of these new variables, 
G{X;k) = G{X-,K + qA), the star product becomes 
e^p{^[K>-dxl - hx^di^l + qhZ'E -vl-qk- V^d^]). In 
the rest of the manuscript, we simplify the notations by 
dropping the tilde. 

For uniform and stationary solutions, the dependence 
on X is lost and the equations reduce to 



D{k) * G"{k) = 1 , 

G^(k) = G^(k) * ^^{k) * G^in)* 



(2) 
(3) 



with D{k) = zu + e{K) — E^(k) and the star product is 



simplified to * = exp {^qld^'E ■ V^— fe • VkS^]). In the 
_E = limit, the star product between Wigner transforms 
reduces to the usual product between Fourier transforms 
and the equilibrium equations are immediately recov- 
ered. The quantum Boltzmann transport theory 16| cor- 
responds to performing a gradient expansion in E. 

Dimensional reduction. When the energy scale asso- 
ciated to the electric field {\q\Ea) is much larger than 
any other relevant energy scale, the dimensionality of 
the system is reduced since the hoppings along the di- 
rection of the field are inhibited as electrons experience 



Bloch oscillations with a vanishing amplitude eo/ 
(see e 



\q\E 



17[ for the non- interacting case). These turn 



into a steady state in the presence of a small scattering. 
The very strong field limit can be worked out analyti- 
cally in Eqs. ^ and ([3]), showing that Green's functions 
(and therefore self-energies as well) lose their dependence 
along the direction of the field (i.e. in Kx) and the prob- 
lem reduces to the equilibrium Hubbard problem in lower 
dimensions [l^ : 

G^{K)^[w + e{K^)-^^{m,K^)y\ (4) 
G^{k) ^\G^{w,K^)\'^^{m,K^_), (5) 

with e{K,±) = (a/27r) J dK^ e{K± + KxUx) the correspond- 
ing lattice dispersion relation. This crossover from an 
out-of-equilibrium problem in d dimensions to an equi- 
librium problem in lower dimensions constitutes one of 
the main results of this work. When the field has sizeable 
projections on several principal axes, those directions are 
similarly suppressed. It extends to any model on a pe- 
riodic lattice under the assumption that a steady state 
can be reached and it is valid prior to any approximation. 
Below, we support this result by following the crossover 
of selected observables in a two-dimensional system. We 
obtain the results by a mean-field approximation for Ej/ 
which we present now. 

Dynamical- Mean Field Theory. — The DMFT is a 
non-perturbative approximation scheme, particularly 
well suited for strongly correlated systems in equilibrium. 
Here, we extend it to our non-equilibrium steady-state 
regime. The self-energy E(7(k) arising from the electron 
interactions is taken to be local {i.e independent of k, 
which is exact in the d — >■ oo limit) and equal to the one 
of an impurity problem associated self-consistently to the 
original lattice problem [l^ . The action of the impurity 
reads in the (a, b e ±)-basis 



^ = EE Utdt'-ci{t)Q-'''\t-t')ci{t') 



J2a dt Uc^^{t)c^t)cl{t)cl{t) 



(6) 



where the impurity non-interacting Green's functions (of- 
ten referred as the Weiss effective fields) are determined 
self-consistently through 



e^(tn) = [G«(tn)-i+E0(n7)] ' 
G^Um) 



e^(t^) = |e^(tn)| 



|G«(.n) 



(7) 

(8) 



with G^/-^'(cc7) = (a/27r)''/dK: G^/-^(k). 
Let us recast Eq. 1^ into 



G' 



G^{k) + G^(k) * AE«(ti7) * G^ 



(9) 



where G^{k) is the non-interacting solution of G^(k) * 
Dq{k) = 1 with Do{k) = m + e{K) — E^ and 
AE^(cc7) = E^(cc7) - E^. E^ is an arbitrary 'dissi- 
pative' constant (ImE^ < 0) which helps the conver- 
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FIG. 2: (color online) (a) Density of states at fixed repulsion {U = 
3) for different values of the electric field. The crosses correspond 
to the equilibrium d = I Hubbard model, (b) Density of states at 
fixed field {E = 3) for different values of the electronic repulsion. 
(T = 0.1,r = 0.5). 
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FIG. 1: (color online) (a) Spectral function (integrated over Ky) at 
(7 = 3 for i? = 3, 5, and 10. (b) Momentum distribution function 
at C/ = 3 for B = 1,3, and 5. (T = 0.1, r = 0.5). 



gence of the integrals below. Once functions are ex- 
pressed in the mixed (r; K)-space (after a Fourier trans- 
form to the real-time domain) one has Gq-{t;k) = 
-iexp [i XT^/adr' e {k + qEr') - irE^j Qir) where e(r) 

is the Heaviside step function. Equations (O and (|8]) 
together with Eqs. © and (O generalize the DMFT to 
our driven steady-state scenario and constitute one of the 
main results of this work. 

Our DMFT algorithm in the steady state goes along 
the following lines, (i) The lattice Green's functions are 
determined. is obtained by solving Eq. ^ self- 

consistently starting from Gq and from a guess 

for the impurity self-energy (typically S^'' := 0). The 
right-hand-side of Eq. ^ is evaluated numerically via a 
discretized form of the following integral representation 
of the star product: 



dr 



' / (r - 

(t';k- 



(10) 



G^ is then obtained similarly from Eq. ([3]). (ii) The 
Weiss fields of the impurity are determined by Eqs. ([7]) 
and ([8]). (iii) The self-energy S(7 is then obtained by solv- 
ing the impurity problem in the steady-state. In princi- 
ple this can be done by any out-of-equilibrium method. 
We simply solve it by means of IPT, treating the on-site 
interaction to the second order in U: 

= -iU6iT)6atg''\T) - U\b I g"''(r) | ^G^' (r) , 

(11) 

This perturbative treatment of the impurity provides a 
good qualitative picture of the many-body physics of the 
lattice problem and has already been used out of equi- 
librium 0. The whole procedure, (i) to (iii), is repeated 
until convergence is achieved. 

Results. — We study the effect of the electric field 
and illustrate the dimensional crossover by implementing 



numerically the previous DMFT equations and by com- 
puting several physical observables as a function of the 
field. We consider the two-dimensional {d ~ 2) Hubbard 
model on a square lattice with nearest-neighbor hop- 
pings: €{k) = eg [cos{Kxa) + cos{Kya)]. We work at half- 
filling and absorb the Hartree shift by expressing Green's 
functions in terms of the variable tu' = tij — U/2. Below, 
we drop the prime. The density of states of the ther- 
mal baths is chosen to be a Gaussian with a width con- 
trolled by W, yielding E^(ci7) = ~iTw{w / ^/ttW) where 
w{x) = [1 — erf(— ix)] e~'^ is the Faddeeva function. E^ 
can be expressed using the FDT (and fio — 0). We con- 
centrate on the effects of the electric field on the metallic 
phase. We work at very small temperature T, small dis- 
sipation r (but large enough to allow for a stable steady 
state), and at fixed hopping amplitude eg. Hereafter, all 
numerical results are obtained with a = q = 1 and all 
energies are given in units of eo. 

Spectral and momentum distribution function. In 
Fig. [ija), we study the out-of-equilibrium spectral func- 
tion p(e, k) = — ImG^(K)|ro=e/7r averaged over Ky. 
This quantity can be accessed experimentally by an- 
gle resolved photoemission spectroscopy (ARPES). In 
Fig. [ijb), we also plot the momentum distribution 
function n(/«) = 1/2 - G^^(r = 0;k). To illustrate 
the crossover from the two-dimensional to the one- 
dimensional system, we show that both observables lose 
their dependence on as the field intensity is increased, 
until they are completely invariant under translations. 

Density of states. For a fixed U, the effect of the elec- 
tric field on the local density of states (DOS) p(e) = 
— ImG'''^(n7 = e)/7r is illustrated in Fig. [2](a). Starting 
from the equilibrium metallic solution, the quasi-particle 
peak is first rounded up for intermediate field intensities. 
For very strong fields, it develops a structure peculiar 
to the equilibrium paramagnetic DMFT solution of the 
d — 1 Hubbard model (the DOS of which is given with 
crosses for comparison). In Fig. [2](b) we follow the lo- 
cal density of states when varying U at fixed E. For 
— U <^ \q\Ea, the double peak structure of the (non- 
interacting) d = 1 Hubbard model is recovered. For in- 
creasing values of U, the onset of the Mott transition 
occurs as the quasi-particle peak decreases and the insu- 
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FIG. 3: (color online) (a) Distribution function for different values 
of the electric field at fixed U = 3. Notice how the distribution 
function starts from a Fermi-Dirac distribution for E = 0, and how 
it ends as a Fermi-Dirac distribution for E = 7 (T = 0.1, P = 0.5). 
(b) Steady-state current density as a function of the electric field 
for different values of the electronic repulsion (T = 0.05, r = 0.25). 



lating lobes grow. 

Distribution function. In equilibrium, the FDT im- 
plies that the distribution function is /^-independent 
and given by the Fermi-Dirac distribution /(e) = [1 + 
gte-z^oj/Tj-i^ p£ equilibrium, it is now /^-dependent 
and we study in Fig. [SJa) the local energy distribution, 
corresponding to the distribution function of the impu- 
rity, </)(e) = i [l-t-G'^(tn)/ImG«(ti7)]^^^. For weak 
fields, electrons located just below the Fermi level are 
excited to states above the Fermi level. This can be seen 
as the system heating up, increasing the heat transfer to 
the reservoirs and balancing the injected power [20| . As 
the field intensity is increased, this pocket of electrons 
sent above the Fermi surface gets broader and the dis- 
tribution function shows two steps at energies of order 
±qEa. These steps get wider but also thinner for a very 
strong field and the distribution approaches the Fermi- 
Dirac one, showing that equilibrium is reached. 

Current. In Fig. ^h) we study the steady-state cur- 
rent density, J cx 2q / dK {n)^ ^.^{k) , as a function 
of the field intensity and U . The establishment of a fi- 
nite current is made possible by two scattering processes. 
One is due to the electronic interaction whereas another 
is due to the coupling to the dissipative bath [see e.g. 
the finite current at J7 = in Fig. [3fb)]. For all values 
of the electronic interaction, the steady-state current is 
first linear with the field intensity with a linear conduc- 
tivity which is a decreasing function of U , as expected by 
linear response theory. After going through a maximum, 
the current is exponentially damped at very strong fields. 
The damping is again consistent with the dimensional re- 
duction occurring at very strong fields: the physics are 
the ones of an equilibrium (current-less) Hubbard model 
in one dimension. A similar non-linear characteristic has 
already been found in the 2d t — J model [2lj where the 
spin degrees of freedom were claimed to act as a dissi- 
pative background. For the case of a Hubbard sample 
strapped between two leads at different chemical poten- 
tials, a maximum in the current- voltage characteristic 
was also found in the case of a two-leg ladder [22] but 
not in the case of a 3d sample 23| , maybe due to a lack 
of a dissipative mechanism in the bulk. 



Conclusion. — We have investigated a correlated 
metal in a non-equilibrium steady-state regime driven 
by a constant electric field. By taking into account the 
need for a heat dissipation mechanism, and by writing 
the Schwinger-Dyson equations in a covariant fashion, 
we revealed a crossover to a lower dimensional system in 
equilibrium as the field intensity is increased. We have 
generalized the equilibrium DMFT equations to the case 
of correlated systems driven by electric field, which allows 
a comprehensive understanding of the steady-state solu- 
tion at a relative low computational cost. By solving the 
impurity in the IPT approximation, we illustrated the 
dimensional crossover on a two-dimensional system. We 
showed that ^-dependent observables lose their depen- 
dence in the direction of the field as this one is increased. 
We repeated the analysis for local quantities and showed 
their convergence to equilibrium. In particular, the dis- 
tribution function was shown to converge towards the 
Fermi-Dirac distribution, completing the proof that it is 
indeed equilibrium which is reached in the d — 1 dimen- 
sions. 
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